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TECHNICAL NOTE No. 1780 

CKARTS FOR THE CONICAL PART OF THE 
DOWNWASH FIELD OF SWEPT WINGS 
AT SUPERSONIC SPEEDS 


By Jack N. Nielsen and Edward W. Perkins 


SUMMARY 

Analytical expressions hare teen derived for the dcwnwash 
throughout the entire induced flow field of lifting triangles of 
infinite chord with leading edges either in front of or Behind the 
Mach cone. These expressions have teen determined from the line- 
pressure source theory of R. T. Jones and the lifting— triangle 
investigation" of H. J. Stewart. Based on these analytical results , 
downwash charts have teen prepared from which the downwash field 
may he determined for any practical combination of leading— edge sweep 
angle and flight Mach number. These charts for the lifting triangle 
of infinite chord are basic to Idle solution for the downwash field 
of any finite swept wing wherein the contributions of the tips and 
the trailing edges are to be determined by the method of lift can- 
cellation as employed by Lagerstrom and others. The conical part 
of the downwash is obtained from the charts, and the nonconical 
part must be obtained by other means such as the lift-oaneellation 
method. 


INTRODUCTION 

One of the principal requirements for a rational analysis of 
the longitudinal stability of aircraft 1s a knowledge of the down- 
wash field behind lifting surfaces. Theoretical methods based on 
lifting— line theory for determining downwash for conventional lifting 
surfaces at subsonic speeds are well known. However, before satis- 
factory agreement between experiment and theory was obtained, the 
theoretical downwash had to be corrected for the local effects of 
the wake and the vertical displacement of the trailing vortex sheet 
(e.g. , reference l). The methods of reference 1 are applicable to 
conventional airplanes throughout the suboritioal speed range by 
use of the Glauert-Prandtl rule (reference 2). No practical method 
exists at present for the calculation of downwash in the supercritical 
speed range. 
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Several methods "based on solutions of the linearized differ- 
ential equation of compressible flow are available for determining 
the theoretical downwash at supersonic speeds. By the use of the 
conical flow theory of Busemann (reference 3 ), Lagerstrom in refer- 
ence 4 has developed analy tical expressions for the downwash field 
of certain uniformly loaded lifting surfaces. By superposition of 
these lifting surfaces, lagerstrom has determined the dcwnwash field 
of a flat rectangular wing and a flat trapezoidal wing. In addition 
he has indicated a method for determ i ning the effect on the dcwnwash 
field of adding a trailing edge to a lifting triangle of Infinite • 
chord to form a f ini te wing. Heaslet and Lomax (reference 5) have 
determined the downwash for points along the intersection of the 
chord plane with the vertical plane of symmetry behind a finite tri- 
angular wing with subsonic leading edges. They have also determined 
an approximate solution for the downwash in the vicinity of this line 
at an infinite distance behind the wing. The Schlicting vortex theory 
■ can be used to calculate the downwash induced by various span load 
distributions. The concepts employed Eire analogous to the Prandtl 
lifting-line theory and as such consider only the spanwise lift 
distribution- and hence neglect any possible effects of chordwise lift 
distribution. 

No determination has been made of the downwash throughout the 
induced- flow field of a swept wing in supersonic flow. The solution 
. for the downwash field for a lifting triangle of infinite chord is 
basic to the solution for the downwash field of any finite swept 
wing wherein the contribution of the tips and the trailing edge are 
to be determined by the method of lift cancellation as employed by 
lagerstrom in reference 4. far a finite swept wing, the conical 
downwash field of the lifting triangle of infinite chord represents 
the entire downwash field except for the regions of influence of 
the tips end trailing edge. Within these regions, it represents 
an appreciable contribution to the downwash. 

The purpose of the present investigation was to determine 
analytical expressions for the downwash throughout the entire induced 
flow field for flat lift ing tr i angles of infinite chord with either 
subsonic or supersonic leading edges. These expressions are used in 
the construction of nondimens ional downwash charts, which, when used 
in c onj unct! on with the lift— cancellation method, provide a practical 
means of determining the downwash field for finite swept wings. Although 
the analysis was carried out for Mq 2 = 2, the charts axe presented 
in a form •mn.lHng them applicable to any supersonic Mach number. The 
charts axe subject to the usual limitations of linesix theory. 
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SYMBOLS 

Primary Symbols 
angle of attack 

real part of a complex function 
free-stream velocity - 

free-stream Mach, rnrnber 

cotangent of the sweep angle of the wing leading edge 

slope of the surface of the wedge section in the 
streamwise direction 


x,y,z 


x'^ySz* 

9 


u,v,w 


W 

w 

t 

B 


e 


B 


O 


e 


longitudinal, lateral, and normal coordinates with the 
origin at the apex and the z = 0 plane coincident 
with the chord plane of the wing. (The positive 
directions are indicated in fig. 1-) 

oblique coordinates, x* = x-my, y* = y-mx, z* = z tJl-uP- 

potential of the perturbation velocities 

perturbation velocities in x— , y— , and z-directions,' 
respectively 

downwash function (w + iw) 
harmonic conjugate of w 
complex variable (se*®) 

y(y/s) g+ (z/g) g 

1+rJl-ij/^) 2 -( z/x) 2 

^ KS 

V l +jy /l-«? 


downwash angle (-w/v o ) 



MCA TH Ho. 1780 


k 


(^) 


E 




couplet e elliptic integral, of the second kind with, 
modulus 

couplex variables 


E(S ,/\/l— e Q 4 ) incomplete ellipti c integral, of the first kind with 
modulus tJl—B'f and sine amplitude £ 


i[\ 


AL 


Vi-t 2 Vi-(i-fio 4 )S 2 J 


E(i , Vl - 6 o 4 ) incomplete ellipt ic integral of the second kind with 
modulus a/i-s 0 4 and sine amplitude | 


C 

T 

5 


r Am 


i-Vh 2 


'•] 


H> Vl -£ 2 

constant of integration 

real part of a complex variable 

imaginary part of a complex variable 

real part of the incomplete elliptic integral of the 
first kind 


imaginary part of .the incomplete elliptic integral of 
the first kind 


E, 


E, 


k 

k* 


sn 

cn 


real part of the incomplete elliptic integral of the 
second kind 

imaginary part of the incomplete elliptic integral of 
the second kind 

modulus of elliptic integrals ( a/i—Bq 4, ) 
comodulus of elliptic integrals ( r\fl— k 2 ) 

i 

Jacobian sine amplitude elliptic function 
Jacobian cosine amplitude elliptic function 
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dn Jacobian delta amplitude elliptic function 

VT sn(F r ,k) 

tjo sn(Fi,k’ ) 

d incomplete elliptic integral of first k-fnfl of sine 

amplitude £ and modulus k 


Subscripts 

P value at point P 

Q value at point Q 

c . value at Macb. cone 

vc value between plane wave and Macb cone 

x refers to conqolex variable £ 

2 referB to complex variable i] 

r refers to real part of a complex variable 

i refers to imaginary part of a complex variable 


LIFTING TRIANGLE OF INFINITE CHORD WITH 
SUPERSONIC LEADING EDGES 

In the theory of supersonic conical flow, the downwash field of 
a flat lifting triangle is symmetrical about both the chord plane «nri 
the vertical plane of symmetry. The downwash field above or below the 
flat lifting triangle at angle of attack a is the same as the down- 
wash field below a nonlifting triangular wing of the same plan form 
with a streamwise wedge-shaped section of half angle a, since the 
flows above and below the wings are independent. (See fig. 1.) The 
downwash field for this nonlifting triangular wing may be determined 
from its streamwise perturbation velocity field which may be obtained 
from the results presented by E. T. Jones in reference 6. The stream- 
wise perturbation velocity field, or the ur-field, for the nonlifting 
triangular wing is the sum of the ur-fields for two line press ur e 
sources coincident with its leading edges. S imi larly, the downwash 
field for the wing is the sum of the downwash fields for the two lire 
pressure sources. However, since the downwash field for one line 
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source is the reflection in the xz— plane of the downwash field for 
the other line source, it is sufficient to determine an expression 
for the downwash of a single line source. The complete downwash field 
is then determined as the sum of the field given by this expression 
and the field given by this same expression with the sign of the 
lateral dimension reversed. 

R. T. Jones in reference 6 gives the following equation for the 
perturbation velocity u at any given point in the field for a line 
pressure source in front of the Mach cone. 



In this equation |dz/dx| is equal to a as previously discussed. 

It should also be pointed out that the positive root of the ra di c al 

I p .0 

+ z should be used in conjunction with the principal value 
of the inverse cosine. In the succeeding analysis, the downwash 
field for the upper surface of the nonlifting wing will be determined. 
However, since the downwash field for this wing is antisymmetric about 
the chord plane, it will be necessary to reverse the sign o'^-ohe 
result for application to the upper surface of the lifting triangle. 

The perturbation potential <p at any point P may be evaluated 
as the line integral of u along the line parallel to the x-axis 
through the point P. 


( 2 ) 


for which the lower limit xq corresponds to a point outside the 
region of influence of the line source where the values of <P and u 
are zero. The region of influence of the line source is bounded by 
the Mach cone and the plane waves from the line source tangent to the 
Mach cone. (See fig. 1.) Figure 2 shows two possible paths of integration 
for equation (2). The path corresponding to point P of this figure 
intersects only the Mach cone; whereas the path correspo nding to point 
Q intersects both the plane wave (y ,2 +z ,2 = 0) and the Mach cone 
(x 2 = y 2 +z. 2 ). Consider first point P: The perturbation potential 

is given by 
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The vertical velocity component is obtained by differentiation 
of cpp . Thus 


dz 


-1: 


chi 

8z 


dx 


w 


z 2 

since at the lower limit u = 0 in accordance with equation (1). 


Tn considering point Q, it is convenient to evaluate equation 
(2) in several steps: first through the plane wave, then from the 

plane wave to the Mach cone, and finally from the Mach cone to the 
point Q. In -pa ssing thr ough the tangent plane, the value of the 
term cob^{x x in equation (l) Junqps from 0 to it; thus, 

the discontinuity in u is finite, and there is no discontinuity 
in cp. Hence, the lower limit of the int egral s may be taken as the 

, z/v/m 2 — 1 

value of x at the plane wave, y/m + — . Between the plane 

wave And the Mach cone the value of u is constant at — V 0 ma/ 1 
so that within this region the potential is 



This equation shows that the vertioal velocity is constant in the 
whole region between the plane wave and the Mach cone. 1 

When the Integration is carried from the Mach cone to point Q, 
it is noted from equation (5) that the potential at the Mach cone is 


Yoina 

90 = ' 




(7) 


and therefore the potential at point Q is 

9„ - - Sg- (jprss - 1 _ i . f z 


VP+z 1 " 


udx 


(8) 
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Differentiating to obtain the vertical velocity 


Wn = 


Q 


= T 0 a ♦ / 




dx 

8z 


(9) 


Equations (4) and (9) give the value of the vertical velocity at 
points P and\ Q. Although both equations contain the same integral, 
the actual value of the integral vill depend on the path of integration. 
The integral is evaluated in Append i x A for each path; thus, for 
point P 


i: 


rj y 2 +z 2 
and for point Q 
x 


Ir 


chi . T 0 a 

dx = cos— ! 

az * 


yy*+z 2 


V^+z 2 Vy’ 2 +z»2 


( 10 ) 




+z‘ 


8u 

— dx = 
8z 


2s“ coS -i ( — zz. 

Tf V 






VqO- (11) 


If the values of the integral given by equations (10) and (11) are 
substituted into equations (4) and(9), respectively, identical results 
for the vertical velocities at points P and Q are obtained. Thus 


T 0 a 

*P = W Q = -IT COB 


,-x( 


yy* + z 2 

Jjz+z 2 Vy , 2 +z 12 


) 


(12) 


Equation (12) represents the vertical velocity due to one line- 
pressure source. The vertical velocity due to both pressure sources 
is obtained by adding to this result the contribution of the other 
pressure source which is determined by substituting — y for y in 
equation (12). In accordance with previous considerations, the 
vertical velocity for the triangular lifting wing is of opposite sign 
to the vertical velocity due to the two pressure sources, thus 


v 0 a r , 

•— — J COS"” 1 

y(y-mx)+z 2 

* l 

. Vy 2 +z 2 -z 2 (m 2 — l) . 


+ cos- 


y(y+mx)+z 2 

■ </y 2 +z 2 V(7+mc)^ (n£— l) 


]} 


(13) 
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LIFTING TRIANGLE OF INFINITE CHORD WITH 
SUBSONIC LEADING EDGES 

- © 

The lift of a flat triangular wing with, leading edges swept behind 
the Mach cone has been treated theoretically by H. J. Stewart in refer- 
ence 7* In that paper a differential equation is derived for the 
complex! downwash function, the real part of which is the vertical 
velocity. To determine the vertical velocity w, it is necessary to 
integrate the differential, equation and determine the real part of 
the solution. 


In the notation of this paper, Stewart’s differential equation 1 b 
dw 2^a (l+£) 2 


mE(Vl-m 2 ) (?Vf /2 (^^ /2 

The actual integration of this equation is accomplished In 
Appendix B, from which it is found that 


W = w+iw = 


gVoOSp 


mE( Vl-m 2 ) (l+s 0 2 r 


[ (l+2s 0 2 ) F(|,Vl-6o 4 ) 


0*) 


+ E(^a/1-So 4 ) +F(n,Vl-So 4 ) -E(q, C (15) 

In this equation, the complex variables £ and rj are related to 
y/x and z/x by the following equations: 


t. = 


(y/x) + 1 (z/x) 

1 +V 1 -(y/ 2 ) 2 -( z A ) 2 


(16) 


E 


j 


(17) 




1 

Vl+So 2 ^ 2 


(18) 


The value of the constant C may be determined from the boundary 
condition that w = —V Q a on the wing. 
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The vertical velocity w is the real part of the downwash 
function W in equation (15). Thus, it is necessary to determine 
the real par£ of the elliptic integrals appearing in the equation. ■ 
The methods used to obtain the real pant of the incomplete elliptic 
integrals of the first and second kinds are described in Appendixes C 
and D, respectively. Using the results of these append ix es, the 
final expression for the vertical velocity is 


w = 


2Y 0 os 0 


mEO/l-mP) (1+Sq 2 )' 


{ 


(l+2s Q 2 ) F^Xi, Vl-So 4 ) 


+ Vl-So 4 ) + 


( 1— So 4 ) ' cri V^i( 1 — ^1 ) 
1— tfi II— (1— s 0 


+ F(^X 2 , V I -8 o 4 ) — E(<v/l 2 ) V l - 6 o 4 ) 

(1— S p 4 ) Gz /Jxz (1— A-g) [1— (l— S q 4 )^] 

1—02 [1— ^ c 1 — 1 

On the basis of conical flow considerations, the vertical velocity 
is constant along rays through the apex of the lifting triangle and 
thus depends only on y/x and z/x. In order to determine the vertical 
velocity for a ray defined by y/x and z/x, the real and imaginary 
parts of 5 and q must be determined in accordance with equations 
(16), (17), and (18). From r x and Si, the real and imaginary parts of 
£, the values of o x anti X x are obtained using equations (C12) and 
(C13). Similarly, values of 02 and \ z are calculated for the 
complex variable q using these same equations. The vertical, velocity 
w for the ray defined by y/x and z/x is then determined from 
equation (19) using these values of 01 , X x , 02 , and X 2 . 



RESULTS AND DISCUSSION 
Downwash Charts 

The equations presented in this report permit deter m i n ation of 
the downwash at any point in the induced flow fields of lifting 
triangles of infinite chord. The equations have been used to determine 
dowpwash c har ts for a n umb er of such surfaces. The range of leading- 
edge sweep angle used - in 'the calc ula tion is sufficient to include 
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all practical configurations, and the number of lifting surfaces for 
whioh the charts hare been determined is sufficient to permit inter- 
polation. The charts are presented ip figures 3(a) to 3(h) and show 
lines of constant de/da within the Mach cone in a plane perpendicular 
to the flight direction. Since the dovnwash field is conical, the 
downwash pattern is similar in all planes perpendicular to the flight 
v direction. Only one quadrant of the Mach oone is shown since the 
downwash pattern is symmetrical with respect to both the horizontal 
and vertical planes of symmetry. 

In determining the charts for the various lifting triangles, 
the values of de/da were calculated at a number of points in the 
Mach cone for Mq = The lines of constant de/da were then 

determined from cross plots. Although the calculations have been 
made for Mq = *[ 2 ., in which case de/d a. depends only on'm, 
y/x, and z/x, it can be shown by the Glauert— P randtl rule that 
de/da for any Mach number depends only on m (y-x/l^ 2 - l)/x, 

and ( Z VM?^L)/x. These parameters have been utilized in the charts 
of figure 3 so that the results are applicable at Mach numbers other 
thany2. , 

In determining the value of de/da at any point in the field 
'for a given finite swept wing, the value of de/da for the 
corresponding lifting triangle of infinite chord is first determined 
from figure 3. If the point in question is not within the region 
influenced by the trailing edge or the tips, this value will represent 
the total de/da. However, if the point lies within the region of 
influence of the trailing edge or tips, the value determined from 
figure 3 represents only part of the actual value of de/da. For a 
finite lifting triangle, the remaining part, which consists of the 
contribution to the downwash of the trailing edge, can be determined 
by the method of lagerstrom (reference h) . In the more general case 
of a finite swept wing, it is necessary to take account of the sweep 
angle of the trailing edge and the effect of the tips by a similar 
method of lift cancellation. 


Variation of Dowhwasb Distribution With 
Leading-Edge Sweep Angle 

The plots of figure 3 show the changes in the downwash distribu- 
tion as the sweep angle of the leading edge is increased with respect 
to the Mach cone. When the leading edge of the wi ng is considerably 
ahead of the Mach oone (large values of m -i/Mq 2 -! ) , ds /da is large 

throughout the Mach oone and the rate of change of de /da with distance 
above or below the chord plane is small for considerable distances from 
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the chord plane. Between the Mach cone and the plane waves from the 
leading edges, which are tangent to the Mach cone, the value of de/da 
is constant and equal to unity. The surfaces of constant downwash 
intersect at the lines of tangency of the plane waves an d the M ach 
cone. As the sweep of the leading edge is increased (m A /M 0 ^—l 
decreases), the line of t ang ency moves along the Mach cone toward 
the chord plane and the rate of chan ge of de/da with distance from 
the chord plane increases. When m^/Mo 2 — 1 approaches 1, the lines 
of tangency and th e leadin g edge approach coincidence on the Mach 
cone; nnfl when m .*/ M q 2 — 1 is- equal to 1, the surfaces of constant 
downwash intersect on this line of coincidence. 

When the leading edge Becomes subsonic, there are two qualitative 
changes in the downwash pattern. First, the surfaces of constant 
downwash intersect at the leading edges rather than at the lines of 
tangency, which no longer exist; second, a region of upwash occurs 
Between the Mach cone and the leading edge. This region of upwash 
increases in extent as the sweep of the leading edge increases. In 
addition, the rate of ohange of de/dtx with distance above or. Below 
the lifting surface increases rapidly, and the downwash near the 
surface of the Mach cone Becomes very small. In the limit, as the 
direction of the 3 aading edge coincides with the free-stream direction, 
the downwash field disappears. 


Additional Factors Affecting, 
the Downwash 

When the charts of this report are used in conjunction with the 
lift— cancellation method to obtain the downwash Behind finite swept 
wings, certain factors which have not Been discussed thus far must 
Be considered. In the linearized theory of supersonic flow the 
assumption is made that the trailing vortex sheet is coincident with 
the extended chord plane. This assumption, although valid for very 
sma3-l angles of attack, wild, probably" cause an appreciable error at 
large angles of attack for which the trailing vortex sheet is dis- 
placed from the extended chord plane due to the action of the downwash. 
An additional effect, which is important in the immediate vicinity of 
the wake Boundaries, is inflow to the wake resulting from the viscous 
properties of the fluid. The influence on the downwash at subsonic 
speeds of these two factors has Been investigated and methods for 
predicting their influence have Been reported in reference 1. However, 
sufficient experimental data are not yet available to permit the 
development of similar methods for Bupersonio flow. 


Another factor that is neglected By the linear theory is the 
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curling up of the Yortex sheet "behind the wing tips. This effect is ■ 
neglected in subsonic flow calculations for conventional high-aspect- 
ratio wings. However, for the low-aspeot— ratio wings which will 
probably be used with some types of supersonic aircraft, this effect 
may be important. 


Ames Aeronautical Laboratory, 

National Advisory Comm i ttee for Aeronautics, 
Moffett Field, Calif. 


appendix: a 

EVALUATION OF INTEGRAL IN EQUATIONS (h) AND (9) 
The integral to be evaluated is 



in which 


u 


— r.p. 


OOBHL \ 

* V*£-i Wy' I 2 +z' 2 / 


(Al) 


(A2) 


The value of chi/dz is determined from equation (A2) ami substituted 
into equation (Al), with t he quantities x’, y’, and z* replaced by 
x-my, y-mx, and z/y/l-m^, respectively. Thus 

I f z tesl ... _ — (A3) 

* “c; 2^2 t (y-mx) 2 -« 2 (ii?-l)] -A 2 -? 2 -* 2 

From the identity 

x-my 1 f y(m 2 — 1 )-z Vm 2 — 1 _ y (m 2 — 1 )+z*J m 2 — 1 

(y—mx) 2 — z 2 (za 2 — l) 2mz .i L y+z /ym 2 -d-'‘~-mx y— z mx 

(A4) 

the integration can be accomplished, by separating the integral into 
two parts . Thus 
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ik 


2 it l 


m(y 2 +z 2 )— i( y+z Vm 2 — I ) 1 

' (y^/m?— 1— z ) V z 2 — y 2 — z 2 


— tan —1 


m(y 2 +z 2 )— z(y— z /y/m 2 — l) 


L(y A /m 2 -l+z) a/ 1 2 — y 2 — z 2 


V 


JX 


(A5) 


7F+^ 


Before substituting the limits, it is convenient to replace the 
inverse— tangent difference by the equivalent inverse cosine. 


V n a I 


r y(y-mx) + z 2 

^"1 

L 1 

_ Vy 2+z2 J x) _ 


}_ 

Vy^+z 2 


(a6) 


In this equation the positive roots of the radicals are to be used 
together with the principal values of the inverse cosine. At the 
lower limit where 


I s = y 2 +z 2 


and 


(x-my) 2 = (y—mx) 2 — z 2 (m 2 -l) 


the principal value of the inverse cosine is zero for y < x/m and 
n for y >x/m. It can be shown from figdre 2 that y<x/m corresponds 
to point P and y>x/m corresponds to point Q. 


Therefore, for point P 

r 


4 y 2 +z 2 
and for point Q 

du 


2s. ta = M OOB V nr' * g — ^ 

dz * V A /^5I^2 a/v ,2 +Z 12 / 


/ 


Vy 2 + z 2 


. "Voa 

— dx = 

dz * 


cos 


yyjJLif ^ 

^Vy 2 +Z 2 Vy ,2 +Z' 2 ' 


-ft 


(A7) 


(a8) 


appendix b 

INTEGRATION OF STEWART ’ S DIFFERENTIAL EQUATION 
In dete rmini ng W from Stewart’s differential equation 
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aw _ gypo. 


2\2 


(i±T) 


^ niE(vT^) (C 2 +s 0 2) 3/2 ^2 lV /2 


B o‘ 


(Bl) 


it is convenient to express the downwash function as the sum of three 
integrals 


W = 


2Y 0 a s 0 £ 


mE( a/ 1— in? ) ( 1 — Sq^ ) 
+ (l_e o 2)2 J . 


1— Sn 4 


I 




L Sq 2 J />/ £^+6q 2 


A at (i-^o 2 ) 2 

r ^ 1 

(r+B0^ S/a Vt 2 =0 2 +l S 0 2 J 

(C’V'+U 3 ' 2 /^ +Bo 5 J 


(B2) 


Introducing the change of independent variable 


5 = 


Vfc 2+s o £ 


the first two integrals hecome 

d i 


f- 

U Al 


I 


*1 


V^ 2 +So 2 V? 2 s 0 2 +l 


= F(|,Vl-eo 4 ) 


(B3) 


(B4) 


' l 


(^ 2 +S 0 2 ) 3 / 2 /^ 2 S 0 2 + l Sq 2 Vl-d- 60 2 )^ 2 

- - , S ° g 4 . ?(S, Vl-^° 4 ) * , 4 T e (S, VW 
(1-60 ) d o 2 ( 1_s o ) 

(B5) 


Making use of the transformation 

1 


n = 


Vi+s 0 2 ^ 2 


(B 6) 
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the -third integral hecomes 

r . . - r 

^ (C 2 s 2 + 1) S/S Vt 2 + s n 2 ^ 


JTfdT) 


//1 - Tl 2 Vl “ TV (1 - So 4 ) 


- |E(T),Vl - So 4 ) - “ Bo 4 ) 

- B~ 4y l L J 

(B7) 


(1 - So 4 ) 


*3 


Substituting these results into equation (B2) the equation for the 
downwash function "becomes 


W = w+iw = 




mE( (l+s c 2 ) 2 


(l+2s 0 2) F(|, Vl-So 4 ) 


+ E( 6, Vl-so 4 ) + r(ri, Vl-a'o 4 ) - E(t), yi-e 0 4 ) 


+ C 


(B 8 ) 


The constant C is to "be chosen so that v, the real part of W, 
is zero on the Mach cone. 


APPENDIX C 

DETERMINATION OF THE REAL PART OF THE INCOMPLETE 
ELLIPTIC INTEGRAL OF THE FIRST KIND 

The Jacohian normal form of the incomplete elliptic integral 
of the first kind in the complex plane is determined as follows: 

P T+1B dg ■ 

F(t + 18, k) ~ / - 7 =-== (Cl) 

J o V 1 -! 2 Vl-^ 2 0 

In determining the real and imaginary parts F r and F i of this 
integral, it is convenient to introduce the Jacobian sine amplitude 
elliptic function. 


sn (F r + iF^,k) = t + 18 


(C2) 
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Malting use of the addition formula for the sine amplitude function 


sn(F +1P k) sn ( :F r> k )cn(tFi,h)dn(lPi,k)-t-sn(iF i ,k)cn(F r ,k)dn(F r ,k) 
r ^ 1— k 2 sn 2 (F r ,k)sn 2 ( iFi,k) 


(C3) 


Substituting the following relationships 




cn(iF i} k) = 


cntF^k’) 






equation ( C3 ) becomes 

sn(F r> k)dn(Fijk t )+i sn(Fi,k’ )cn(Fi,k’ )cn(F r ,k)dn(F r ,k) 


sn(FjvflFi,k) = 


cn 2 (F^,k* )+k2sn 2 (F r ,k)sn 2 (F ; i,k t ) 


Separating equation (C5) into its real and imaginary parts 


t = 


sn(F r ,k)dn(F^^k* ) 


cn^F-pk* j+k^n^Fyjk^n^F^k’ ) 


sn(Fpk* )cn(F^^k t )cn(F r ,k)dn(F r ,k) 
' cnSfFpk* )+k 2 sn 2 (F r ,k)sn 2 (F i ,k* ) 


(C5) 


(c6) 


(C7) 


Introducing the substitutions 

X = sn 2 (F r ,k) 

a = sn 2 (Fi,k T ) 


(C8) 
(C 9) 


and making use of the identities 
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cn(F r ,k) = \fl-\ 
dn(F r ,k) = 

cn(Fi,k* ) = \J±-<3 

dn(Fi,k» ) = v/l— k ,2 a 
equations (C 6 ) and (C7) become 

T 2 = X(l— k , 2 g) 

(l-g+k^g ) 2 

52 = g(l-g) ( 1 -^) (l-^X) 
( 1— cr+k^cr) 2 


(CIO) 

(Cll) 


Simultaneous solution of these equations yields 


= [ (1+t 2 +5 2 ) — \/(1+t 2 +6 2 ) 2 - 4t 2 ] {l+k 2 ^ 2 -^ 2 ) - 4l^r g } 


lH^r 2 


(t 2 +6 2 -^) 


(t 2 +5 2 -X) - [Xk^T^+e*) - 1] 


2/_2.c2> 


(C12) 

(C13) 


Therefore, to find the real and imaginary parts F r and F^, of 
the incomplete elliptic integral pf which the sine amplitude is r+i5, 
it is first necessary to determine X and g from equations (Cl 2) and 
(C13). Then from equations (C 8 ) and (C9)> it follows that 


(C14) 


F r = F(v/X,k) 

F ± = F(Va,k«) 


(C15) 
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APPENDIX D 


DETERMINATION OF THE REAL PART OF THE INCOMPLETE 
ELLIPTIC INTEGRAL OF THE SECOND £IND 

% 

The Jacobian normal form of the incomplete elliptic integral 
of the second THnfl in the complex plane is given as follows: 


Tt-i8 ■ . — 

E(T+i8,k) = f ~rr- d£ 
J o jn 5 


(Dl) 


In deter mining the real and imaginary parts Ep and E.^ of 
this integral, it is convenient to introduce the transformation 

I = sn(p,k) = snCpp+ip^k) (D2) 


where the real and Imaginary parts of p Eire p and p^, respec- 
tively. With this transformat ion equation (Dl) becomes 


Fr+IFi 


E(r+i6,k) = / dn 2 (p,k)dp 

J o 

where the upper limit 

T+i§ 

FxH-lFi = f 


■k 2 ! 2 


(D3) 

m 


is evaluated in accordance with the methods of Appendix C. 

In integrating equation (D3) it is convenient to perform the 
integration in two steps. The first integration is along the real 
axis to F r , and the second integration is aiong a line parallel to 
the imaginary axis from p = F r to p = F^iF^ Thus equation (D3) 
becomes 


F r I 

E( T+i5,k) = J' dn 2 (p r ,k)dp r +iy^ 


n 


dn 2 ( F r +ip i , k ) dp i 


(D5) 


\ 



In determining the real part E? of E(r+i8,k) from equation (D5), It Is noted that the first 
Integral Is real and is equal to E(snF r ,k). 


8 


Tne determination of the real part of the second integral of equation (Dp) may he performed 
with the help of the addition formula for the delta amplitude elliptic function 


dn(F r +i4J.£,k:) 


dn(F r ,k)dn( lpj,k) — k^snCFj jk) s a( jFi > k)on(E r ^k)cn( lFj,k) 


i i_ 'i O / jm 1_\ 


(D 6) 


Using the relationships for sn(iFi,k), cn^iFijk) and dn(iFi,k) (See Appendix C.) in this 
equation, and substituting the imaginary peart of the resulting expression for dn(F r +iF^,k) into 
the second integral yields 


r.p, i 


i . 

dn 2 (F r +J 4 i.j_,k)dp^ = 2k 2 Bn(F r ,k)cn(F r ,k)dn(F r ,k) 


1 cn(p 1 ,k , )dn(p 1 ,k t )sn(p 1 ',k T ) 
[on 2 (qi,k' )+k s Bn 2 (F r ,k)Bn 2 (p i ,k’ ) ] 2 

(DT) 


By direct integration of this equation 


np. 


i 


i 

dn 2 (F r +ip i/ k)du i 


k?sn(F r , k) cn(F r ,k)dn(F r , k) 


sn 2 faA') 

1— Bn 2 (Fi,k’ )dn 2 (F r Vk) 


(D8) 


n VJ I . -L.T n ..14. n jy J.V A j)4 un 4. a n 

uuiiiuj.nl life ouo roouj-UD iur oun j.jlj.oo emu oouvuu. uiia^icuo 


Er 


, ^ 2 aa(?r^)cn(F r ,k)dn(F r ,k)8n 2 (Fi,k») 
r> 1— sn 2 (Fi,k* Jdi^tFr^k) 


(D9) 


With the substitution of equations (C8) and (C9) of Appendix C, equation (D9) becomes 
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E r = E(VT,k) 


kVl a/I-1 Vl-k 2 ! a 
l-oCi-^x) 


(DIO) 


Thus the real part of the incomplete elliptic integral of the 
second kind, of sine amplitude r+iS is obtained by evaluating X 
and a from equations (C12) and (C13) and substituting these values 
into equation (DIO). 
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Figure t. — Wing and Mach wave system 
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Figure 3 . — Continued 
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Figure 3. — Continued 
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Mach line 
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/V — Mach line 
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^ ^ M§ - / tan cu * .4 


Figure 3 . — Continued 
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